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Abstract: We embed a holographic model of an U (1) charged fluid with Galilean invari- 
ance in string theory and calculate its specific heat capacity and Prandtl number. Such 
theories are generated by a i?-symmetry twist along a null direction of a A/" = 1 super- 
conformal theory. We study the hydrodynamic properties of such systems employing ideas 
from the fluid-gravity correspondence. 
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1. Introduction 

Systems at strong coupling continue to be one of the least tamed areas of modern physics. 
This is particularly puzzling considering the number of physical situations that seem to 
be described in terms of such models. Examples include non-BCS superconductors at the 
quantum critical point, the hydrodynamics of strongly coupled microscopic theories or the 
quark-gluon plasma (QGP). More curious still is the fact that until relatively recently we 
have even lacked a description of strongly coupled toy models, for instance, planar A/" = 4 
SYM with a large 't Hooft coupling. 

Notably, hydrodynamics has been the subject of intense study and, while there have 
been significant achievements in this area such as those of Navier, Stokes and Kolmogorov, 
there remain many phenomena, like turbulence, that still lack a full theoretical description. 
Importantly, an ability to completely describe the generic behaviour of fluids would lead to 
a deeper understanding of a vast range of models. This is because a typical feature of field 
theories with a long-wavelength expansion is a sector well described by the hydrodynamic 
regime. 

For a fluid description whose microscopic origin is strongly coupled we cannot apply the 
usual perturbative methods to calculate approximate transport coefficients. As such, these 
fluids represent an interesting theoretical challenge. Moreover, beyond purely theoretical 
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considerations, there are several practical applications where controllable models of such 
fluids would be useful. One example concerns the dynamics of the QGP whose transport 
coefficients are now accessible to experiment [1, 2, 3]. The interest in applying the gauge- 
gravity duality to this system comes from the fact that the correspondence describes a 
large class of ideal fluids and the QGP seems to be approximately ideal [4]. Using the 
general lessons learned from the application of AdS/CFT has already led to a qualitative 
improvement in understanding the properties of this strange state of matter like its low 
viscosity to entropy ratio [5]. 

The AdS/CFT correspondence [6, 7, 8] provides a useful tool via which we can learn 
general lessons about certain strongly-coupled systems from their gravity duals. To date 
the most extensive use of this conjecture has been in the study of non-abelian gauge theories 
since early examples of the correspondence related A/" = 4 SYM, and similar models, to 
Einstein gravity. Like several dualities it relates a theory with a large expansion parameter 
to another theory with a small perturbation parameter. However, unusually for such a 
correspondence, the "master field" corresponding to the small curvature regime lives in 
one higher spacetime dimension than the strongly coupled theory, indicating a holographic 
description. 

This paper is principally concerned with the fluid-gravity sector of the AdS/CFT 
correspondence where it was discovered that the restriction of certain large N gauge theories 
to a long-wavelength regime is dual to a simplified gravitational description in which the 
relevant quantities are determined completely by hydrodynamic conservation laws [9] . This 
programme of studying duals to the fluid description of strongly coupled conformal field 
theories began with the seminal works of [10, 11] where calculations of the graviton retarded 
Green's functions were made at the linearised level of gravity. A significant achievement in 
[12] was a procedure to extend the previous calculations perturbatively to higher orders in 
derivatives of fiuid velocity. It should be noted that the results obtained by these methods 
can give only qualitative predictions about nature because, as yet, no observed phenomenon 
is known to have as an underlying description a large N gauge field. 

Recently, significant effort has been expended in studying the non-relativistic limit 
of strongly-coupled field theories and their gravity duals. The interest in these theories 
in part lies in the fact that it has proven difficult to find generalisations of the AdS/CFT 
correspondence to other asymptotics. One of the few cases where this has been man- 
aged are the "asymptotically Schrodinger spacetimes" where the dual field theory satisfies 
the Schrodinger algebra [13, 14, 15] possibly with a central extension. In the light-cone 
coordinate system the bulk metric of the ground-state of these spaces has the form: 



where /3 is some number. For the case above /3 can be scaled out by a boost but it has 
non-trivial effects outside of the ground state. Based on the suggestions of [16, 17] it later 
proved possible to embed these spacetimes in string theory [18, 19, 20] although there 
as yet remains issues of interpretation of boundary quantities [21]. For other works on 
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Schrodinger spactimes and non-relativistic CFTs see [22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 
32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49]. 

It seems reasonable to ask whether we can extend the fluid-gravity correspondence 
to discuss fluids with Schrodinger symmetry and perhaps shed some light on outstanding 
problems in non-relativistic fluid mechanics. The case of an uncharged fluid has already 
been considered in [50]. However, adding charge to a fluid has previously led to new con- 
ceptual insights such as the necessity of a parity violation term for a relativistic fluid, with 
a gravitational dual description, in a charged background [51]. Our focus is then naturally 
drawn towards understanding charged fluids with Schrodinger symmetry. In recent papers 
[52, 53] an C/(l)-charged, asymptotically Schrodinger spacetime at finite temperature has 
been discussed and some linear transport coefficients calculated. To obtain these results a 
truncation of the 10-dimensional effective Type IIB string action was isolated. Specifically 
we begin with a stack of D3-branes rotating in [54]. The Kaluza-Klein reduction of 
this solution is associated with a RN — AdSs black hole. By applying a Null Melvin Twist 
[55, 56], or alternately a TsT transformation [57], to the 10-dimensional metric we induce 
Schrodinger symmetry on the boundary spacetime. Here we would like to extend the anal- 
ysis of these charged, thermal systems to hydrodynamics at first order calculating all the 
relevant transport coefficients. 

A fluid with Schrodinger symmetry whose gravitational dual was {d + 3)-dimensional 
would occupy d spatial dimensions. For comparison, an alternate approach to achieving 
non-relativistic symmetry is given by the "Galilean conformal algebra" which produces 
fluids moving in d+1 spatial dimensions. Intuitively this corresponds to permitting super- 
luminal communication and is achieved mathematically by taking a suitable scaling of 
hydrodynamic variables and spacetime coordinates. This approach preserves the form of 
the equation of state from the relativistic theory to the non-relativistic theory but projects 
out sound waves whose dispersion relations take the form oj (x k giving an incompressible 
fluid [58]. Unlike the Schrodinger algebra it does not admit a central extension; this can 
be interpreted as meaning the model describes gapless non-relativistic fluids. We do not 
consider this approach here because, as was discussed in [58], the charge conservation 
equation to first order reduces to demanding incompressibility. Hence the charged fluid is 
a trivial generalisation of the uncharged fluid under this type of scaling. 

Our aim is to detail the first order corrections to a charged fluid with Schrodinger sym- 
metry. We demonstrate how the general form of these corrections can be obtained by 
reducing conformal, relativistic currents to their non-relativistic counterparts. We then 
consider a five-dimensional, asymptotically Schrodinger, charged black hole spacetime and 
construct its fluid dual. The corrections up to first order in boundary derivatives of velocity 
to the metric, dilaton, gauge field and massive vector field are calculated. Finally, using the 
holographic dictionary in Schrodinger spacetimes, we shall calculate the boundary values of 
the stress tensor and gauge field, obtained from our charged black brane, to determine the 
dependence on charge and temperature profiles of the non-relativistic transport coefficients. 

The paper is organised as follows. In the second section we shall discuss the theoretical 
formulation of non-relativistic hydrodynamics. This discussion will not be pedagogical but 
instead highlight a few important specifics necessary for the development of the rest of the 
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paper, in particular, the structure of the hydrodynamic stress tensor and charge current at 
first order in derivatives. In the third section a map is constructed that hnks the transport 
coefficients of the relativistic charge current to their non-relativistic counterparts. This is 
an extension of previous work, notably, we seek to preserve the mappings of [50] between 
relativistic stress-energy-momentum (SEM) tensor transport coefficients and their non- 
relativistic counterparts. In the fourth section we shall construct an action principle for 
the charged, asymptotically Schrodinger black hole and determine its field content. For 
this calculation we closely follow the work of [52]. We then, in section five, briefly describe 
the method of determining the transport coefficients by solving the Einstein equations 
order by order in derivatives of velocity as developed in the seminal paper [12]. Again the 
discussion of this section will avoid many of the heavier details and instead highlight the 
most relevant features. We refer the interested reader to [12] for the minutiae. Having 
obtained the asymptotically AdSs metric and gauge field to first order we take the TsT 
of this solution to find the corresponding metric and fields with Schrodinger symmetry. 
Finally we calculate the boundary values of these fields and determine the non-relativistic 
transport coefficients. 

2. Hydrodynamics 

In this section we shall formulate the basics of non-relativistic hydrodynamics. Our starting 
point is to consider the fluid as a thermodynamic system. We then indicate how transport 
coefficients are related to the derivative expansions of the spatial stress tensor, energy 
current and charge current. In particular we discuss conventions, such as which operators 
to keep, in the effective field theory and set out future notation. 

We shall assume our fluid is in local thermodynamic equilibrium in the neighbourhood 
of any point. The defining thermodynamic potential is the Gibbs potential: 

G = G{P{x),T{x),N{x),Q\x)) (2.1) 

where P{x), T(x), N{x) and Q\x) are the local pressure, temperature, particle number 
and charge of the fiuid. For an isobaric and isothermal process in an open system it 
represents work done as particles move from one region of local equilibrium to another. 
Changes between equilibrium states of the system are specified by: 

--(i) -Hi) "^Hi) "^Hf) '''^''^•^> 

\'^^/T,N,Qi \'^^ / P,N,Qi X*^^^ / P,T,Qi V'^V / T,P,N 

^ V ' ^ V ' ^ V ' ^ V ' 

V{x) -S{x) IJ,{x) lJ-q,l{x) 

where V{x), S{x), and Hgj^x) are the volume, entropy and chemical potentials respec- 
tively. It will be more profitable to consider volume densities of thermodynamic quantities 
which sets G{x), S{x), N{x) and Q\x) to g{x)^ s{x), n{x) and q\x). We shall also 
consider only a single U{1) charge, the generalisations being clear, and hence can drop 
the charge group index /. Upon integration (2.2) yields an equation of state satisfied by 
the system whenever it is in equilibrium. Importantly, in a transition between equilibria. 
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the components of dG complete the specification of the thermodynamic state of the fluid. 
However, from a macroscopic perspective, the intensive quantities are a priori unknown. 

The thermodynamic coefficients do not indicate how quantities flow between different 
patches of local equilibrium. This is the realm of non-relativistic hydrodynamics^ and as 
such we need to supplement our knowledge of the extensive variables with the local fluid 
velocity v^{x) and mass density p{x). The hydrodynamic regime is characterised by four 
conservation equations: 



d+p + di {pv') = 

d+ {pv') + djU'^ = 



d+qnr + diJl 








(2.3) 
(2.4) 

(2.5) 

(2.6) 



where we have used + to denote the time coordinate to match our later interpretation of x"*" 
under light-cone reduction. These equations are the continuity, momentum conservation, 
energy conservation and charge conservation equations respectively of the fluid. To zeroth 
order in derivatives of velocity and temperature we can expand the undetermined tensor 
objects, ILij, Jl and j'^^ as: 



n'^' = pv'v^ + P^^6'^ 

fe = (enr + Pnr + ^pv'^ I V" 



Jnr 



QnrV 



(2.7) 
(2.8) 
(2.9) 



where enr and Pnr are the fluid's energy density and pressure. We shall call (2.7), (2.8) and 
(2.9) the stress tensor, energy density current and charge density current at zeroth order. 
These quantities characterise a perfect fluid which does not lose energy due to internal 
friction as it has no viscosity. 

Following [9] we consider adding terms to our spatial stress tensor and currents with 
single derivatives of velocity. Our undetermined tensor quantities in (2.4) and (2.5) take 
the form: 



1 . 

Je = ( e„r + Pnr + Ij^pV 



-wS'^^dj In 



<rp fl ^lq- 

' T' T . 



) ri{T,J^,t^) 



3 At ' + 



{^2 (rp a _i_if3 



(rp jl /iqN 

: T' T ) 



(2.10) 



(2.11) 



^Several quantities such as relativistic and non-relativistic pressure are traditionally labeled with the 
same symbol P. We shall have cause to investigate both relativistic quantities and their non-relativistic 
counterparts. To avoid confusion, unless an object is explicitly defined as non-relativistic, for example 
the mass density p, we shall adopt a notation where the subscript nr indicates non-relativistic while an 
unlabeled object is relativistic. 
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where we have assumed a flat background. The stress tensor Iljj is unchanged from the 
uncharged case of [50] but the energy current j* has received an additional correction which 
vanishes when the charge is set to zero. In these expressions we have used the following 
quantities: 



r+ix,y,z) 



X \ 2 

y. 



1 + A 1 + 



32 z2 



(2.12) 



assuming our non-relativistic fluid occupies d spatial dimensions. The coefficients of (2.10) 
and all but the last coefficient of (2.11) have standard physical interpretations; r]nr repre- 
sents the fluid's shear viscosity, C„r the bulk viscosity and the thermal conductivity [59]. 
We shall call the new coefficient, ro, the contribution to the energy current from charge. 

As regards determining the first order corrections to the current vector we note that 
at zeroth order the conservation equations (2.3)-(2.6) can be written as: 

9+p + v'dip + pdiv' = 

d+v^ + v^dy + -d'Pnr = 
P 

d+enr + di {enrV') + PnrdjV^ = 
d+qnr + v'diQnr + gnrdiV' = 

where the equation of state for the fluid relates €nr and Pnr- Hence, if we obtain a complete 
solution to the above equations they allow us to replace time derivatives of the variables 
€nr, p, Qnr and for spatial derivatives at first order making an error in our final results 
that, overall, is second order in derivatives and can therefore be ignored. We can thus write 
our non-relativistic current vector as: 



Jr. 



(2.13) 



Here Knr is the non-relativistic diffusion constant and I3nr tbe parity violation coefficient. 
We shall call the tensor objects 7nr and Fnr the contributions of energy density and mass 
density to the charge current respectively. 



3. Light-cone reduction of charged relativistic fluids 

One way of obtaining the form of the first order corrections to (2.7), (2.8) and (2.9) is to 
light-cone reduce the SEM tensor and charge current of a relativistic fluid. In particular we 
shall consider a conformal, relativistic fluid and then the light-cone reduction will lead to a 
hydrodynamic system with Schrodinger symmetry [50] . We begin this section by stating the 
conformally invariant, relativistic conserved currents, specifically, our choice of variables. 
We then summarise the previous work of [50] on the uncharged fiuid before determining a 
map between relativistic charge coefficients and their non-relativistic counterparts. 
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Determining the relativistic conserved currents provides the starting point for our 
analysis. Using conformahty, the choice of Landau frame and relativistic invariance, as 
shown in [9], it is possible to decompose these currents to first order in derivatives in the 
following manner: 

T^"" = {e + P) u^v" + Pr?^'^ - 2?] (^P^"P'^^V(«n^) - ■^^^V«n°P^'^^ (3.1) 

f = qu^ - KqP^'^Vuq - iP^^V^e - Uea/i-y ''u'^V^u^ (3.2) 

where P^'^ = u^u'^ + rj^'^ . The decomposition in (3.2) is not explicitly conformal in Kq and 
7 and, therefore, they must satisfy the following relation: 

{d+l)Kgq + {d + 2)je = Q (3.3) 

The final term of (3.2) is specific to four dimensions as this is the only case where there is 
a first order parity violating contribution. 

The parity violating term of (3.2) may at first be alarming. Traditional approaches 
to determining the corrections to such as Israel-Stewart theory [60, 61] typically ignore 
or argue this coefficient away. However, as was shown in [51], in a theory with a charged 
background, consideration of triangle anomalies leads to its presence. By demanding posi- 
tivity of the entropy current it is possible to find an expression for 13 in terms of the gauge 
anomaly coefficients and the hydrodynamic variables. We refer the reader to [51] for the 
general case and shall be content here to determine 13 for the specific case of a fluid dual 
to an asymptotically Schrodinger, ?7(l)-charged black hole. We shall find it to be fixed by 
the charge q and the fiuid pressure P up to an anomaly dependent constant. The other 
terms in (3.2) have more sedate interpretations as the charge diffusion coefficients Kq and 
contributions of the energy density to the charge current 7. 

As was demonstrated in [50] at first-order in derivatives of fluid velocity there ex- 
ists a map between the relativistic SEM tensor variables, {u^,e, P,r]), and non-relativistic 
(i)*, p, Cnr, Pnr, Vnr) Variables. We shall seek to maintain these relations and augment them 
with our charge variable maps. In particular we can use them in our reduction of the 
charge current. Summarising the results of [50] we begin by assuming that our fluid lives 
on a Minkowskian background with metric: 

ds'^ = 2dx+dx^ + dx^ (3.4) 

and that the relativistic hydrodynamic variables and velocities depend only trivially on the 
x~ direction. This suggests we make the following identiflcations: 



e-nr + ]^pv^ ] (3.5) 





= p 


rp+i 






-{ 


rp — i 


= -Ji 
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which come from comparing SEM tensor conservation equations in our choice of coordinates 
and (2.3)-(2.5). The component of the SEM tensor implies the following identification 
between variables: 



p = (e + P) (n^ 
while the T"*"* component indicates that: 



(3.6) 



V = — I diu^ 

'+ p 



u 



diP 



P \ (e + P) 

The other quantities of importance we list for completeness: 





= pv' 


Pnr 


= P 


Vnr 


= rju 








-V 


Jl 





nr + Pnr + l^pV 



"^Vnr Pnr 
P 



3dj€nr djp 



2 tr 



P 



(3.7) 
(3.8) 
(3.9) 

(3.10) 
(3.11) 



where (3.7) has no bulk viscosity term as conformal invariance of the parent relativistic 
theory set this to zero. For a conformal relativistic fluid in {d + 2) spacetime dimensions the 
equation of state is supplied by tracelessness of the SEM tensor and implies e = (d + 1) P 
and therefore, using the above maps, the non-relativistic fluid satisfies enr = ^Pnr- 

In [50] the final two terms of (3.11) are eliminated in favour of diT and the resultant 
coefficient is interpreted as the thermal conductivity kt- This can be done because the 
equation of state for an uncharged fluid is given by: 



/i 



d+2 
2 



where a is a constant as detailed in [34] . However for a charged fluid there is an additional 
scale in the problem, Hq, and therefore our equation of state takes the more generic form: 
Pnr = T^g (^, ^) which prevents us from eliminating dienr and dip completely. We 
shall return to the interpretation of these terms in section 5. 

Considering now the charge density current it is clear that it satisfies the conservation 
equation: 

+ dd' = 

Using the expression for jf^ from (3.2) and the maps (3.6)-(3.11) we find that has the 
form: 



J 



qu+ + U{u+fd^v''ejk 



(3.12) 



where we have used the scaling relation (3.3) to annihilate the term proportional to 6 and 
-€ij. It is a satisfying occurrence that the only correction to the identification 



set ej 
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j"*" = at first order is a piece wiiich accounts for the anomalies in the relativistic theory. 
Indeed, in the holographic model we shall construct in future sections, if the Chern-Simon's 
coupling in our action is set to zero, then qnr is just a scaling by u+ of q. Reducing the 
spatial part of the current leads to: 



3 



J V 



d + 2 



+ 



+ 



3 
P 



{d + 4) r/n+ 



d p 



u~ 
djp 
2p 



6'^ + Ue'^ 



2e 



+ 



r\,q 



Additionally, in light of the above expansion, it seems reasonable to define 
and I5nr = 'U{u^)'^. The remaining terms, proportional to diP and dip, merit further 
consideration for the same reasons as the final terms of (3.11) and we shall return to them 
in light of section 5. However, using the equation of state to find e in terms of e^r the 
spatial part of the current can be rewritten as: 



d 



e'^v^idjVk-dkV,) + v'e^''djVk 



d + 2 



+ 



Qnr 
1]nr 



{d + 4) r]„ 
d c 



5'^ + 



+ 



P 



sii 2 (d + 2) I3nr^nr 

■' H , e 

a p 



2 {d + 2) Unre 
d p 

2p 



nr ^ij 



'2c n; 



+ 



which represents the maximal possible reduction into our chosen non-relativistic operators. 
This matches (2.13) if we identify the coefficient of djCnr with ^nr and that of djp with 

r nr- 

So far we have only related the non-relativistic and relativistic variables. A priori we 
have no knowledge of the functional form of the relativistic charge coefficients in terms of 
local charge and temperature profiles. Hence we also do not know the functional forms 
of the non-relativistic coefficients in terms of the corresponding non-relativistic charge 
and temperature. In the next section we specialise our discussion to one particular fluid; 
a charged, Schrodinger invariant fluid which has as its dual a charged black hole with 
asymptotically Schrodinger boundary. The holographic dictionary in Schrodinger space- 
times will then readily allow us to calculate these dependencies; an otherwise difficult 
problem for a general fluid. 



4. Charged, asymptotically Schrodinger black hole construction 

As our focus is now on Schrodinger fluids with gravity duals we can apply the holo- 
graphic dictionary in Schrodinger space-times [18, 19, 62] to our investigation. This pro- 
vides a map from an asymptotically Schrodinger space-time to a boundary field theory with 
Schrodinger symmetry and hence allows us to determine expectation values of the spatial 
stress tensor, energy current and charge current at zeroth order. In particular, we begin 
with an uplift for RN — AdSs which we can TsT transform to give the bulk fields boundary 
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Schrodinger symmetry. To calculate the temperature and charge profiles corresponding to 
our dual black hole we will need to construct a suitable action yielding these bulk fields. 
We lean heavily on the formalism of [52] and find a 10-dimensional action from an effective 
description of Type IIB string theory whose Kaluza-Klein reduction leads to a five dimen- 
sional theory with the correct field content. It turns out that this reduced action contains 
only four fundamental fields; the gravitational field, the dilaton, an i?-charged gauge field 
and a massive vector that is now standard fare in Schrodinger spacetimes [16, 18, 19, 20]. 
We complete this section by calculating the thermodynamics of our bulk spacetime, in 
particular, the specific heat at constant pressure, particle number and charge. 

We shall take m to be the mass of our black hole and Q to be its charge. From [63] a 
suitable RN — AdSs ansatz is: 

-r'^f{m,Q,r)dt^ 



+ # + A(i 



(1) 



+ 
2 



r^f{m, Q, r] 

2 

Ao 1 + cRl 



+ [dx^ + dy^ + dz^) + 



^(5) 


= 2(l + *io 




= — 




2 


Fq 


= dAQ 


Q,r) 


m 





(4.1) 

2 \ 1 " 

(4.3) 
(4.4) 

, . (4-5) 

where *5 and are the five- and ten-dimensional Hodge stars defined in the appendix 
and the unusual factors of — are present so that our Aq matches that of [64] . Note that 
by choosing Aq = a{r)dt we only have an electric field turned on. The final two terms 
of the ten-dimensional metric are the five-dimensional metric on the unit given by a 
U{1) fibration over CP^ with a gravi-photon a{r)dt turned on. The one-form ^(x) ^-nd the 
two form J^2) are called the Kahler potential and form respectively and are defined by the 
following equations: 



2 

- J(2) A J(2) 



The fields of (4.1)-(4.5) satisfy the equations of motion from the following truncation of 
Type IIB string theory: 



S 



1 



10 



167rG 



- j (*iol)e-'* + 4(a<I>)2) - 1^(5) A *ioF(5) 



(4.6) 



where we have chosen the string frame. As usual, self-duality of the RR 5-form must be 
imposed after variation to complete the specification of the equations of motion. 

Using the five-form Bianchi identify di^(5) = we find our choice for the decomposition 
of the RR 5-form implies that Fq satisfies the following equation of motion: 



d *5 Fq - ^KcsFq a Fq = 
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where ncs = ~2v^ Chern-Simon's coupling of the gauge field Aq. This follows from 

an action of the form: 



which, upon compactification, we shall use to substitute for the F(5) terms of (4.6). There 
is one caveat in this replacement and it is that the cosmological constant in the five- 
dimensional theory receives contributions from the volume form pieces of (4.2). The correct 
cosmological constant can be determined by examining the equations of motion. 

Also, we now generalise to arbitrary kqs which will clarify the role of the Chern- 
Simon's term in generating the parity violating coefficient 15 of (3.2). Compactifying we 
can set the resultant scalar field associated with the metric, cr, to be equal to the dilaton 
Further, we find $ = is a consistent solution to the equations of motion and hence 
our action (4.6) reduces to: 



^5 = / voIa/ {r^''^ + 12) - 2Fq a hFq + ^^Aq a Fq a F(^ 



(4.7) 



in the Einstein frame where G5 = — t7%t. This matches the Einstein-Maxwell action used 
to calculate the fluid metric up to second order in [64]. 

We would now like to induce Schrodinger symmetry on the boundary of the RN — AdSs 
spacetime given by (4.1)-(4.5). To do this we can apply one of a pair of solution generating 
techniques at the level of the equations of motion; the NuU-Melvin twist [55, 56] or TsT 
[57]. They both take an RN — AdSs x x' manifold where is Sasaki-Einstein and yield 
an asymptotically Schrodinger, charged black hole with a deformed x^- The former, NMT, 
begins by boosting our solution along one of the spatial isometry directions, say y by a 
rapidity of 7. Two T-dualisations along y are then performed with a twist of the one-form 
dtp —7- vdip sandwiched between them. We then boost the resultant fields by —7. Finally v 
is scaled to zero and 7—7-00 while keeping (3 = \ve"' constant. The latter technique, TsT, 
involves a twist in the x~ = i^{y — t) direction by adip between two T-dualisations along 
the ip direction. This second technique can be applied to any spacetime with a C/(l) x U{\) 
isometry. Moreover, when one of the U{1) isometrics is null it can be shown that the two 
techniques coincide [50]. 

Applying a TsT along the Hopf direction ij: as discussed in the appendix, where we 
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also give the heavier details of our notation, leads to the following fields: 



idslo)" 



^(3) 
^(5) 



exp(2$'' 



Q, r){dt + dyf - f{m, Q, r)dt^ + dy^ + fcdx^ 



dr2 



r'^f(m, Q, r) 
~k 



+ 



dip + ^(1) - -^Aq 



{dy + /(m, Q, r)dt) A ( + ^(i) - ^ A 



^(5)+5(2)AF('^) 
1 



(4.8) 
(4.9) 

(4.10) 
(4.11) 

(4.12) 



where we have set q = 1 so that boost and twist parameters of NMT and TsT respectively 
coincide. We have defined: 



k = l + /3V(l -/(m,Q,r)) 



(4.13) 



in correspondence with the notation of [18] and note that our results are in agreement with 
[52]. The operation *"q is the Hodge dual in the Melvinised spacetime and its definition is 
also detailed in the appendix. 

The self-duality of the 5-form F"^^ allows us to determine a relation between *5Fq 
in the un-Melvinised spacetime and the quantities /, Am and Fq in the Melvinised 
spacetime^. Hence we can write: 



^5 



1 



levrGr, 



J volM»e-2* + 16 - 4e2*) - 4e*F A *'^F - ^e'^^FM A *'^Fm- 



4e-^AM A 4Am 



{Am a Fq) a 4 {Am A Fq) 



4e 
2 



"* ^ ^Fq + F A Am] a 4 (^Fq + FAA 



V3 



M 



-e'FQ A 4'Fq + ^Aq a Fq a Fq 



(4.14) 



where, as was discovered in [52], the equation of motion for the F field is completely 
algebraic: 



F = e 



-2<J> ,// 



1 



Am a *'5' — Fq + F a Am 



and it is therefore only an auxiliary field which is merely present to simplify our action. 



^See appendices for further details. 
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{dsly = [-(3\^f{m, Q, r){dt + dyf - f{m, Q, r)df + dy^ + kdx^] + 



Considering (4.14) we can see that the Kaluza-Klein reduction of the fields (4.8)-(4.12) 
in the string frame is: 

~k 

^r^fim,Q,r) ^^'^^^ 
Aq = ^9.dt (4.16) 

Am = ^{dy + f{m, Q, r)dt) (4.17) 

where the massless one form and non-trivial dilaton field come from the metric while the 
massive vector field originates in the NS-NS two form B^2)- We note that the existence of 
a charged, massless one-form whose boundary value shall be interpreted as sourcing the 
charge of our fluid and a massive vector field which lacks a corresponding conserved current 
on the boundary. 

With the full 5-dimensional zeroth order metric (4.15) available we can now determine 
the thermodynamics of our fluid. Given m and Q, our metric has a horizon whenever 
/(m, (5, r) = and hence will have an associated Hawking temperature T. We take r_|_ to 
be the location of the outermost horizon which turns out to be Killing as our spacetime is 
stationary. The temperature can be found by determining the surface gravity, k, from the 
following formula: 



1 



(4.18) 

r=r_|_ 

where is the null generator associated with the horizon. We are working in the Einstein 
frame attained from (4.15) by conformally rescaling the metric with the dilaton: 



■2, 



ds% = {dsif (4.19) 

For (4.19) the (Killing) vector generating the horizon is proportional to {dt)"' . To determine 
the constant of proportionality we note that in lightcone coordinates, x"*^ and (9+)"" is 
the generator of boundary time translations and hence if we set its coefficient to be unit 
we find: 

r = ^ {dtT 

= (d+r - 2^ (d.r (4.20) 
Substituting this result into (4.18) the temperature, T = is: 

In the limit that Q — this coincides with the result of [50]. Of note, the temperature of the 
relativistic precursor theory (4.1)-(4.5) is given in terms of the non-relativistic temperature 
(4.21) by PT. 
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The entropy associated with (4.15) can be calculated using the Hawking formula: 

s 



4Gb 



where A is the area of the event horizon at r = r+ . The solution turns out to be independent 
of P which is to do with the fact that our metric was generated by a series of boosts and 
dualities from (4.1) as discussed in [56]. For our solution, taking V3 to be volume of the 
horizon, we find the entropy in the boundary is given by: 



4Gf 



V3 (4.22) 



whose form in terms of r_|_ is unchanged from the uncharged case discussed in [50] . 

In (4.20) we have a generator of time translation on the boundary (d^)^ and an 
additional generator ((?_-)'^. The former corresponds to the Hamiltonian H in the dual field 
theory while the latter should be interpreted as the particle number generator N. Hence 
the coefficient of ((?-)" in the null generator is the particle number chemical potential: 



1 

'2/32 



However, this is not the only chemical potential in the thermodynamics of our solution as 
the charge can also vary. This charge chemical potential can be found from the asymptotic 
values of the "boundary time" component of the gauge field and is given by: 

fiq = A+(r+) -^+(00) 

~ 4/3r2 

V6Q(-At)^ 



(4.23) 



Note that we have explicitly defined both the chemical potentials as non-relativistic quan- 
tities. 

We now have sufficient information to specify our density matrix and can therefore 
obtain the thermodynamic potential of our ensemble. The existence of two chemical po- 
tentials implies that the density matrix has the following form: 



exp 



H — fld^ — flqJ^ 

T 



The trace of the density matrix gives us the partition function whence the thermodynamic 
potential is determined by: 

G{V2, T, fi, ^xq) = -T In H(T, ^, ^Xq) (4.24) 

where H = tr (p) and V2 is the two-dimensional spatial volume of the Schrodinger theory. 
This free energy represents the work done by the system when the chemical potential 
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difference between two neighbouring regions of equilibrium changes in an isothermal and 
isochoric process. The Gibbs potential (2.1) can be obtained from (4.24) by a Legendre 
transform. Before we do this however we need to obtain the charge and particle number 
from G. 

We can obtain the relativistic energy E by computing the ADM mass of (4.8). Using 
the relativistic equation of state e = 3P and (3.8) gives us the non-relativistic pressure Pnr- 
Finally, using the integrated form of the first law of thermodynamics and the definition of 
the free energy^ implies: 

G = E -TS - fiN - figj+ 

= -PnriT,fl,fiq)V2 

IGvrGs 1^ + 3 /i +y 

where Ax~ has been introduced to characterise the period of the compactified a;~-direction 
[50]. Using the thermodynamic relations supplied by: 

V5^/t,m,m, \dTJv,f,,^^ \d^^Jv,T,^, \^^^qJv,T,^^ 

V ' ^ V ' ^ V ' ' V ' 

-Pnr{x) -S{x) -N{x) -J+{x) 

we find the following additional quantities: 

jV = _2P^ j+ = ^m/^x^V2 5 = g|Ax-U2 

In particular the entropy, at constant t, matches that obtained from (4.22). Neither 
nor fiq attain their relativistic values on setting /3 = 1 however the combination iJLqJ~^ does. 
We also note that the energy density calculated by inverting (4.25) and dividing by the 
spatial volume V2 satisfies the equation of state for a non-relativistic fluid: 

' " levrGs \ + 3 + J 

— p 

— ^ nr 

which acts check of our thermodynamics. 

The Gibbs potential is given by the Legendre transform: 

G(P„,, T, N, J+)=G + PnrV2 +I^N + fiqj+ 
= flN + flqj+ 

'^Our formula for the free energy does not appear to agree with the on-shell action of [52]. An important 
consistency check of the on-shell action is that one recover the correct entropy i.e.; S = — + l) ^. It 

is not clear to us that the expression in [52] passes this check. On the contrary, because we are working in 
the grand canonical ensemble and have derived (4.26) by calculating the pressure from the ADM mass of 
the black hole, G definitely passes this test. 



(4.25) 
(4.26) 



where all the quantities on the right-hand side are functions of P, T, N and J'^ . It is easy 
to check we get (2.2). Using (4.26) we find: 



^ _ Ax-V2 



4 , 16Atg 2 



3 fi 



SttG, 

The specific heat at constant pressure, particle number and charge is then given by: 

TT^T 



'+ I-' ' T' r y 



2 nA ,A "I ~^ 



4 /i^ 64 

^ ^ n II \ I 



(4.27) 



(4.28) 



3 /ir^ (r, y , J? ) 9 /i^r^ (T, , 5? ^ 

where r+ (T, ^, ^) is given by (2.12). In the limit that the charge goes to zero this 
approaches the result for the heat capacity given by [50]. We note that thermodynamic 
quantities like (4.27) and (4.28) will also apply to the first-order Schrodinger fluid of the 
next section. 

5. A non-relativistic fluid with Schrodinger symmetry 

In this section we go beyond the holographic dictionary in Schrodinger space-times as dis- 
cussed in the previous section and obtain first order corrections to the non-relativistic fluid 
with Schrodinger symmetry discussed above. We begin by applying the fluid-gravity pro- 
cedure of [12] to an asymptotically AdSs charged black hole with translationally invariant 
horizon as done in [64, 65]. We then indicate how the same method could be applied to our 
asymptotically Schrodinger black hole metric (4.15). However we are saved from actually 
carrying out the latter procedure as we can simply TsT transform the, already calculated, 
asymptotically AdSs, charged black brane with first order fluid corrections. Once this has 
been accomplished we calculate the boundary stress tensor complex and charge current 
and identify expressions for the coefficients of (2.10), (2.11) and (2.13). 

5.1 Review of the charged, relativistic fluid-gravity correspondence 

We shall begin by reviewing the fluid-gravity correspondence for a relativistic fluid. The 
reason for doing this is because the procedure for relativistic fluids is more elegant than its 
non-relativistic counterpart making it easier to become aquainted with the ideas involved. 
First, we note that there are several solutions to the equations of motion of an Einstein- 
Maxwell theory with negative cosmological constant that are asymptotically AdSs. A 
solution of interest to us is the "charged, boosted black brane" which, written in terms of 
in-going Eddington-Finkelstein coordinates, has the following form: 

dsl = -2u^dx^'dr - r^f (m^°\Q^^\r^ Ui,u^dx''dx'' + r'^P^ydx^dx" (5.1) 
2r2 

where we have defined: 



Aq = ^^^^u^dx'' (5.2) 



dx^ = , i-dr + Vidx"-) 
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and introduced ^'^^ to emphasise constant objects. This metric has the nice feature that it 
is everywhere non-singular except at r = 0. Note that r plays the role of boundary time 
and as such the boundary metric has the form: rj = diag(— 1, 1, 1, 1). 

Solution (5.1), (5.2) in and of itself is not of interest for this paper as it corresponds to 
a stationary (modulo Lorentz boosts) fluid. Consider instead promoting the constants m, 
Q and to "slowly-varying" fields of the boundary coordinates m{x°'), Q{x") and v^{x°'). 
By slowly varying we mean the following; consider expanding these fields about x = 0. In 
a local patch about this point whose size is set by the scale of thermodynamic fiuctuations 
from local equilibrium, denoted L, we can use coordinate redefinitions and velocity field 
gauge choices to set m(0) = m^^\ (0) = and Q(0) = Q^. It can then be seen that 
higher derivative corrections to ■m{x'^), Q{x°') and Vi{x") are suppressed by higher powers 
of 7^ where T{x°^) is the local temperature of the fluid. We should also note that the 
non-zero temperature and velocities have broken the 5'0(4, 2) symmetry in the boundary 
theory and thus we can loosely identify these fields as Goldstone bosons^. 

The line element and gauge field with promoted constants has the form: 



Ac 



-2n^(x")dx^dr - r^f{m{x''), Q(a;°), r)u/,(x")u;.(x")dx^dx^ + 
+r'^P^^{x'')dx>'dx'' 

2r2 



2 -u^{x")dx'' 



(5.3) 
(5.4) 



Generically the corresponding metric does not satisfy our equations of motion. However, in 
[64] terms of higher order in were added to the above ansatz and treated as perturbations 
upon it. It was discovered that if the m{x"), Q{x°') and v^{x°') fields additionally obey the 
following constraints: 

5.m(x")L=o = -^m(x") d,v\x%^^^ (5.5) 
a,m(x")L=o = -4m(x-) drV,{x^)\^^^ (5.6) 
5rQ(x")L=o = ^^v\x%^^^ (5.7) 



then the Einstein-Maxwell equations, given by variation of (4.7), are satisfied at first order. 
This procedure can be extended to higher orders but our concern here lies only with first 
order corrections. 

To detail this process of solving the Einstein-Maxwell equations more carefully we 
note that the metric (5.3) and gauge field (5.4) can be expanded about the point = 



*The fluid stress tensor wiU be expanded in terms of these Goldstone fields so it can be partitioned into 
collections of terms with the same factors of ^kr ■ 
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to give: 



2dTdr -r'f[ m(°) , r ) {drf + r^dxidx' - 



+ 



-Ix^'d^vfUx'dr - 2x^'d^vf\'^ 



(0) n(o) 



1-/ m^^^Q 



idr) 



dx'dT + 



Ac 



V3 
2 



(5.8) 
(5.9) 



where as previously mentioned we have used coordinate scaHngs and a suitable gauge choice 
for the velocity field to set: 



u{x^') 



-dr + x^d^v'fUx' + . . . 



Note that ultra-locality of our field equations allows us to extend our solution away from 
x'^ = to the entire manifold [12]. Our metric (5.1) and gauge field (5.2) have spatial 
5*0(3) symmetry and therefore we can parameterise the corrections to (5.3) and (5.4) in 
the following form: 

^(1) = -3h^^\r)dTdr + '^—^ {drf + 3r^h^^\r) [dx^f + 



+2r^ 



Q 



r 

(0) o(o) 



j^\r)dTdx^ + r'^a\j\r)dx^dx-' 
dx' 



V3Q(°) .(1), , (1) 

2^2 3i y^)-9i (0 



(5.10) 
(5.11) 



where the first three terms of (5.10) are scalar, the fourth term vector and the fifth term 
a symmetric 2-tensor with respect to SO{2>). Similarly the first term in A^q is a scalar 
and the second term a vector. If we denote the Einstein equations by Eab and the Maxwell 
equations Ma then when Eab and Ma are evaluated on the metric and gauge fields with 
first order corrections they can also be decomposed with respect to (5*0(3). We can further 
separate the Einstein-Maxwell equations in each of these sectors into two groups; constraint 
equations^, which are obtained by contracting the tensors Eab and Ma with the vector 
dual to the one-form (dr)"", and dynamical equations. The constraint equations depend 
on the particular nature of the fluid we are considering in that, as they correspond to 
covariant conservation of SEM tensor and charge currents (5.5)-(5.7), we need to solve 
the hydrodynamic equations for temperature and charge profiles when given a velocity 
distribution. However, the dynamical equations can be completely solved by choosing the 
correct dependence of the undetermined functions in (5.10) and (5.11) on m, Q and r. 



^To call the projected equations constraints is a slight abuse of terminology since we are not dealing 
with an initial value problem here. 
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A relativistic and gauge independent expression for the charged black brane with AdSs 
boundary found in [64] to first order in derivatives is: 



ds\ = —2undx^dr — r'^f{m, Q, r)u^Updx^dx'^ + r^P^^^dx^dx'^ — 
-2rn^ (^u^VxUv^ dx^dx" + {v xu^^ u^u^dx^dx" + 

+2 — Fi(m,Q,r)an^dx^dx'^ ^ Unl^dx^dx^ — 

2 

-12Q^F2(m,Q,r)n^ (p^^Va + Su^Va^^) Qdx^'dx" (5.12) 



\/3Q 
2r 

d 



mr^ ^ 2r^ 



^^F2{m,Q,r) 



P;Va + S-u^VAn^j Qdx^ (5.13) 



where the first line of each definition is zeroth order in fluid derivatives, the subsequent 
lines first order and: 



/(m,g,r+) = 
Fi{m,Q,r) = dx 



X [x" + X + l) 



^ (X + 1) + X2 - ^ 



i^2(m,g,r) = - 1-- + 5, / dx- ^ 1+ + 



These results were obtained by perturbatively solving the equations of motion from (4.7) 
with trivial dilaton field using the fluid derivative expansion method of [12]. It has been 
shown that this solution can readily be uplifted to a 10-dimensional solution of Type IIB 
by adding a suitably deformed term to the metric and folding an Aq term into the 
standard expression for the Ramond-Ramond five form [63]. 

5.2 Charged non-relativistic fluids 

With the well-established relativistic fluid-gravity derivative expansion as a skeleton we 
can now consider moving beyond zeroth order for a fluid with Schrodinger symmetry. In 
particular we shall find the derivative expansion corrections to (4.15) at first order. We 
begin by noting that (4.15) is not regular across the future horizon. One of the nice 
properties of (5.3) was its regularity at all points other than r = 0. This can be remedied 
by translating the x^ and x~ coordinates in the following manner: 

dx~^ — )■ dx~^ o ^ ^ — ^dr 

r'^ j[m, Q, r) 

^ ^ 2/3r2/(L,Q,r)^^ 
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Our metric, gauge field and massive vector field correspondingly become: 



rV {dx+Y + 13^1- f{m,Q,r)) {dx-f + 



1 - f{m,Q,r) 
4/32 



+ {1 + f{m,Q,r))dx^dx + ( ^ + 2/3 ) dx+dr 



2/3 



dx dr + 



Am 



2 r2 

/3r2 
~k~ 



,dr'^ 



2/3 



(3dx 



dx^ 



(1 + /(m, Q, r)) ^ + (1 - fim, Q, r)) ^dx 



dr 



(5.14) 
(5.15) 
(5.16) 



where we have used a gauge choice to remove a dr term from Aq. It should be noted that 
because of the dominant scaling of the x~^ term in (5.14) we lack a well-defined boundary 
metric for the asymptotically Schrodinger, charged black hole. 

Our boundary theory has Galilean symmetry and by boosting (5.14) we obtain a class 
of solutions with the same thermodynamics but non-zero velocity. The boosted solutions 
will be classified by five constants /3, u*, Q and m however these contain no new physics. 
Instead, to obtain results other than the ideal fluid, we need to promote the parameters m, 
Q, /3 and to functions of the boundary coordinates and find first order fluid corrections 
to our fields. A local Galilean boost has the form: 



X 



X -|- v(x)x 



X — )■ X + v(x) • X -|- -V {x)x 



Consider the vicinity of a point x'^ = and use global Galilean invariance to set the 
velocity at this point to zero. Performing our local boost on (5.14)-(5.16), to one derivative 
in velocity, we find the following additional terms: 



.(1) 



2/32r 



2^2 



A 



(1) 



4f = +^(l-/(m,g,r))x.dv 

Again, generically, the metric, gauge field and massive vector field with the above terms 
will not satisfy the equations of motion from (4.14) and as such we need to find suitable 
corrections. Our precursor fields (4.15)-(4.17) have 5*0(2) spatial invariance and thus we 
can parameterise our corrections with respect to this symmetry much as we did with the 
5*0(3) symmetry in the relativistic case (5.10), (5.11). However we face two additional 
complications for the asymptotically Schrodinger spacetime. Firstly, we not only have new 
equations of motion for the massive vector field and dilaton to satisfy but our equations 
of motion have become more complex. Secondly, the 5*0(2) symmetry is not as helpful 



k 

2 r2 

/32r2 



(1 — /(m, Q, r))dx x -|- -(1 + /(m, Q, r))(ix~''x + — -^xdr + 



/3r2 



kx^ 



-dx 



/3x • dv 
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in the Schrodinger case as the 5*0 (3) was in the relativistic case. For example, in the 
Schrodinger case there are two possible vector sectors in the metric coming from dx~^dx^ 
and dx~dx^ compared to one, drdx^, for a relativistic fluid. 

It is clear then that solving the equations of motion from (4.14) to first order in deriva- 
tives would be a cumbersome process. Instead, as previously mentioned, we can perform a 
TsT transformation of (5.12) and (5.13) to obtain the U{1) charged, Schrodinger fluid at 
first order. As we previously assumed trivial x~ dependence in our hydrodynamic variables 
in order to light-cone reduce, we can drop all x~ dependence in the metric coefficients. This 
ties in nicely with the fact that to use the TsT solution generating technique detailed in 
the appendix we require x~ and to be isometry directions of the metric. In particular 
notice that x~ is a null isometry direction in the boundary theory and so the TsT with a 
twist along this direction will coincide with an NMT of our fields. 

Using the identities from [50], in particular the zeroth order current conservation 
equations, our results indicate that the dilaton now has the form: 



= k 



= 1 + f3\^ [1 - f{m, Q, r)] + ^^Z^^' Kcse'^d.v, (5.17) 

where i,j run over the spatial directions {x, z} and we have taken u'^ = /?. The full metric 
at first order is given by: 



(dsl)" = -2u^dx^'dr - r'^f{m, Q, r)u^u^dx^'dx'' + r^P^^dx'^dx" - 

-2ru^ (^u^Vxu^^ dxf'dx" + (Vxu^^ u^.Uydx^dx'' + 2^Fi(?n, Q, r)a^ydx^dx'' - 

_^V^!^csCf ,^^f,^^u _ i2Q!-F2(m,Q,r)ua ( P^V x + 3u^V xu^] Qdx'^dx" - 
mr* V / 

-k{AM)^{AM),dxV (5.18) 
The massive vector field, in the light-cone coordinate system, has the form: 
1 



(1 - /(m,Q,r)) + ^KcsQ^/S^ejfc^^'^^M dx^ + {r^dx+ - pdr) + 



+ ( /3r2 (1 - /(m,Q,r)) + ^^^^f g^.^gi^M u^dx^ + 



+ 2 Fia^a 7 la - 6/3(5^^2 VaQ - -Q T^^" 

+ ^ ) dx^ 

2 £nr ~l~ Pnr J 

where a G and Aq is unchanged from (5.13) by the TsT. Note that all but the 

last line of the metric (5.18) occurs in (5.12) so that all the deformation comes from the 
vector field Am- This makes sense in light of the fact that the massive vector field is the 
only additional dynamical field between the Melvinised and un-Melvinised solutions. 
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5.3 Transport coefficients of non-relativistic fluids 

Now that the bulk metric and vector fields are known we can consider the fluid side of the 
fluid-gravity correspondence. To calculate the hydrodynamic and charge coefficients we 
need to determine the boundary value of the conserved currents associated with the bulk 
metric and the gauge field. However, we face a particular problem in Schrodinger space- 
times due to the slow asymptotic fall-off of the modes. We follow [19], [21] and [50] and 
interpret the SEM tensor of the asymptotically AdSs theory prior to the TsT transforma- 
tion as a tensor complex (collection of fields) in the non-relativistic theory. This complex 
is given exactly by the identifications in (3.5). 

The boundary SEM tensor and charge current for the RN — AdSs black hole have the 
forms: 

T^u = P {r]^u + 4:U^Uu) - 2r]af_,u 
f = qui" - Kg {P>"'V^q + 3gn^V^n^) - Ul^" (5.19) 

The coefficients, in terms of G5 which can be related to the central charge of the field 
theory, are: 

e = 3P P = V = 

1 ~ 47rG5 '^1- \ imr+ ) ' " ^ - 2P ^ 

where we have applied the SEM conservation equation at zeroth order in derivatives: 

n^'VuU^ + = 

e + P 

to the penultimate term of (5.19) to write it in terms of e. Note that our parity violation 
coefficient 13 is indeed determined by the Chern-Simon's parameter kcs^ the charge q and 
the fluid pressure P. 

To convert our relativistic results into their non-relativistic counterparts we need to 
fix the normalisation of our boundary velocity n^. Fortunately we have already isolated a 
suitable choice in (4.20) when we fixed the coefficient of ((?+)" in the horizon null generator, 
1^", to be unit. Hence we take: 

«+ = /3 

= Pv^ - ^S^' (a,/3 - ^-A_5^.p^ (5.20) 

Using the maps (3.7)-(3.11), (3.12) and (3.13) it is now possible to determine all the non- 
relativistic quantities in terms of m, Q, f3 and w*. Modulo the subtlety involving the thermal 
conductivity which we shall discuss next the zeroth order coefficients are: 



Cnr — Pnr Pnr — IGttG"- P 



2P„ 



\/3/3Q 
47rG,, 



1 H ejk&'v'' 

m 



Ax~ (5.21) 
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The apparent disparity for qnr with the thermodynamic result is due to the previous de- 
composition of the combination fiqJ~^ into /ig and which shuffled a factor of two between 
them. At first order we also have: 



inr 



' nr 



Ax" 



Knr — (^4ml3r+) 



2P„r ^i^r 



(5.22) 



1 



4e^ 



1 

2p 



, Vnr 
) — 

P 

inr \ l^nr ~l~ 



Knr ]d'^ +A 
Tjnr 



^nr^nr 
1 

P 



.13 



P 



§ij + ^ '^nr-enr ^ 
P 



Ax" 
Ax" 



where the Ax~ factors were introduced to ensure the above quantities are volume densities 
with respect to the two-dimensional spatial volume V2. The parity violating term I3nr 
already multiplies an object that is first order in fiuid derivatives, see (3.13). Hence when 
expressing it in terms of qnr we have dropped any additional velocity derivatives. Similarly 
for replacing 7 in -fnr- 

We would now like to extract the thermal conductivity. Re-expressing the final two 
terms of (3.11) as the differential of a logarithm and using (4.21) and (5.21) we can write: 



3 

. '^VnrPnr ^ijg^^^ I P,?r 
P \ P 



^Vnr Pnr 

PT 

2,7] jir 



In 1 + 



21n 1 



2ri 



(5.23) 



If we interpret the thermal conductivity as the coefficient of the term with no explicit 
charge dependence then kt has the same functional dependence on Tjnr, Pnr, p and T as 
in the uncharged case of [50] : 

^Vnr Pnr 



pT 



(5.24) 



As promised there is a new term which vanishes if the local charge density is set to zero. 
On substituting for Q in terms of jiq using (4.23) and rearranging we obtain the final term 
of (2.11) with w = Kt- 

We would now like to calculate the Prandtl number for the fiuid. We first note that 
the kinematic viscosity, which compares the importance of viscous to inertial forces, is 
defined by: 

V 



Vnr 
P 



where p is representative of inertia. The thermal diffusivity is defined by: 



X 



PCp„r,N,J+ 



where kt measures heat flow from a region of local equilibrium while pcp^^ j+ measures 
the ability of the region to adjust its temperature to match its surroundings. Thus when x 
is large the region in question quickly responds to the temperature of neighbouring regions 
and equilibriates its temperature. 
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Figure 1: A graph of the Prandtl number against ^iq and T. We note that the Prandtl 
number clearly tends to one when the charge vanishes. 



The Prandtl number is given by the ratio of kinematic viscosity to thermal diffusivity: 

Pr = - 

X 

and thus represents the relative importance of viscous effects and heat conduction in reach- 
ing steady state flow. Using (4.28), (5.21), (5.22) and (5.24) we find this number to be: 

-1 



Pr 



ur^ (t ^ ^ 
+ \ T' T 
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4 



where r_|_ (T, ^, ^) is given by (2.12). Disappointingly this indicates that the fluid does 
not achieve a universal value in the presence of a conserved electric charge unlike the 
uncharged case of [50] where it is identically one. Of note is the fact that Pr is independent 
of the particle number chemical potential and compactification radius Ax". Figure (1) is 
a diagram giving an indication of the dependence of Pr on charge chemical potential and 
temperature. 



6. Discussion 



We began our study by laying the foundations for non-relativistic hydrodynamics in rea- 
sonably general terms. We then demonstrated how to obtain the hydrodynamic derivative 
expansions of a fluid with Schrodinger invariance at first order from a parent charged, con- 
formal, relativistic theory. Specifically we have shown how to generalise the maps of [50] 
to the case of an U{1) charged fiuid at first order. 

We then focused upon the hydrodynamic limit of a particular three dimensional, non- 
relativistic conformal field theory and its dual solution which is an asymptotically Schrs, 
charged black hole. Using the TsT technique on an RN — AdSs precursor we constructed 
an action whose equations of motion had the desired black hole as a solution and isolated 
the thermodynamics. We then obtained first order corrections to the corresponding black 
brane by using a method based on the derivative procedure detailed in [12]. Because of 
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its relative simplicity we first showed how to find corrections to our initial RN — AdSs 
ansatz before discussing how the same method could be applied in the asymptotically 
Schrodinger case. Although in principle we could then have computed the first order 
corrections to our asymptotically Schrodinger, charged black hole using the fluid-derivative 
procedure we noted that it would be particularly cumbersome to do so and hence opted 
instead to use the TsT technique. Thus we arrived at expressions for the metric, gauge 
field, massive vector field and dilaton to first order in derivatives. 

With the corrected fields to hand in principle it was possible to calculate the asymp- 
totic values of the metric and gauge field directly to determine their corresponding con- 
served boundary currents. While this is simple to do for the gauge field ambiguities in 
asymptotic fall off of the metric necessitated that we interpret the boundary SEM tensor 
in the precursor asymptotically AdSs theory as a tensor complex of the Schrodinger invari- 
ant theory [19] . With these identifications it was relatively simple to apply the holographic 
dictionary in Schrodinger space-times to compute the boundary coefficients and with a 
little work obtain the Prandtl number. An important result discovered here was that the 
universal value of one for the Prandtl number of an uncharged fluid no longer holds when 
there is an additional non-zero charge. This suggests that it may be interesting to under- 
stand the consequences of a scaling where the non-relativistic charge and particle number 
were related as this would naturally be interpreted as the charge being carried by the fluid 
particles. We leave this for future work. 

Although our study concentrated on a fluid occupying two spatial dimensions, the 
derivative expansions (2.10), (2.11) and (2.13) apply in any dimension with the caveat that 
the relativistic one-derivative parity violating term only exists in four dimensions. Similarly 
the generalisations to multiple U{1) charges or indeed different internal symmetries seems 
clear and we can determine the hydrodynamic coefficients if we can find a suitable dual black 
hole with the required asymptotics. A more significant limitation of this paper is refiected 
in the fact that our fluid is both conformal and incompressible. To remedy this situation 
we would need to show that we can decouple density fluctuations. As was mentioned 
previously under a suitable scaling it is possible to reproduce the incompressible Navier- 
Stokes equations [58] by suppressing sound modes. However this limit has an entirely 
different symmetry group and the charged and uncharged solutions are trivially related. 
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A. Mathematical conventions 
A.l Hodge duals 

Following [52] we note that the ten-dimensional Hodge dual on (RN — AdSs) x can be 
restricted to the 5-dimensional RN — AdSs manifold, 1 dimensional flbration coordinate 
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and CP^ in the following manner: 

*10 = (_l)(5-'^5)n4+(5-n5)ni+(l-ni)n4 ^^^^ 

where ns, rii and 714 are the number of indices in each part. In particular: 

*io(l) = ^e*volRN-AdS5 A (^-0 + -4(1)) A J(2) A J(2) 
*5(1) = VolRN-AdSs 



*i(l) =e* (^dV + AD-^^g) 



*4(1) = 2*^(2) ^ J{2) 
*4.{J(2)) = J{2) 

where V0I55 = *i(l) A *4(1) when <I> = 0. 

After Melvinisation the Hodge dual on the asymptotically Schrodinger, charged black 
brane spacetime, denoted RN — Schrs , is not equal to that on RN — AdSs and we must 
determine its effect upon our volume forms and gauge fields. It can be shown that objects 
whose terms all contain dx~ pick up a factor of e* when acted on by the Melvinised Hodge 
dual. After Melvinisation the following important objects Hodge dualise in the manner 
shown: 

VolRN-AdSs — 6~*VolRN-Schr5 



hFq = -2e-* 4 (^-^Fq + F a Am) 



The transformation of *5Fq was determined by considering the fact that B'^^^) ^ -^(3) is 
precisely the quantity that needs to be added to make F(5) self-dual with respect to the 
Melvinised metric, see [52]. 

A. 2 TsT transformation 

First note that generically we can write the 10-dimensional metric and five-form as: 

dsfo = g {dx~Y + 2g_adx~dx°' + Qapdx^'dx'^ + (d^ + ^(i) + Aqf + dT,l 

F(5) = A {dx^ A yl(3) + 5(4)) + dx~ A C(4) + -D(5) 

where a, (3 belong to {r, +,2;, z}. As the TsT only ever performs algebraic operations on 
the "0 and x~ isometry directions we only need to keep track of these terms. 

We shall need to T-dualise our solution twice so it makes sense to define a standard 
form for the relevant fields (as in [52]). In particular, we isolate all the dependence on ip 
in our fields and write them in the following manner: 



= 94,^ {d^J + g^^))^ + ... 
B(2) = (^dip + ^g(^)^ A S(^) + . . . 
F(p) = [dtp + 5r(^)) A F(p)^ + F(p)^ 
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where /p indicates the piece of the field with no tp components. We also choose to denote 
the dilaton by <I>o- The T-dualisation of these objects then yields: 

) ^ (-^w) + ••• 



{dsloY = 


1 






= 









2 



We shall denote T-dualised quantities with a '. 

As was stated in the main body of the paper the TsT transformation is formed from 
the following sequence of operations: 

1. T-dualise along the ip direction, 

2. twist along x~ sending it to x~ + aip where a is a constant, 

3. and finally T-dualise along the ■0 direction. 
Applying these operations to the fields of (A.l) we obtain: 



(dlp + An^-^An] ^2 
{dsl,)" = dsl + ^ ^\ ^ ^ + rfSl - ^ {g__dx^ + g.^dx^) 



k ■ " k 



B'l^^ = I {g__dx- + g.a.dx^) A {diP + A(^i) - ^^q) 



^(3) = "^(3) 
^(7) = ^(2) A ^(5) 



where k = 1 + Q?g^^ and ds\ is the original five dimensional metric. From the above 
formulae we can readily identify Am to be jg^udx^. 
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